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The way tension propagates along a chain is a key to govern many of anomalous dynamics in 
macromolecular systems. After introducing the weak and the strong force regimes of the tension 
propagation, we focus on the latter, in which the dynamical fluctuations of a segment in a long 
polymer during its stretching process is investigated. We show that the response, i.e., average 
drift, is anomalous, which is characterized by the nonlinear memory kernel, and its relation to the 
fluctuation is nontrivial. These features are discussed on the basis of the generalized Langevin 
equation, in which the role of the temporal change in spring constant due to the stress hardening is 
pinpointed. We carried out the molecular dynamics simulation, which supports our theory. 

PACS numbers: 36.20.Ey,87.15.H-,83.50.-v 


I. INTRODUCTION 

The fractional Brownian motion (fBm) is a class of 
stochastic process to describe the anomalous diffusion [l| . 
Being a natural extension of normal Brownian motion, it 
can be expressed in terms of its incremental sequence 

x{t) = [ ds r]{s) (1) 

Jo 

with a:(0) = 0, and the fractional Gaussian noise r]{t) 
with zero mean and the correlation 

( 2 ) 

hence, the mean square displacement (MSD) {x^{t)) ^ 
where (• • •) represents the statistical average. 

Many examples of fBm can be found in systems in 
cells and various soft matter systems. A list of Exam¬ 
ples includes the motion of small particle (colloid^ lipid 
granules, etc.) in cells and polymer networks @,i, the 
lateral diffusion inside lipid membranes ii , the pro¬ 
cess of polymer translocation through a nano-pore i 
DNA/RNA hairpin formation 0 and the telomere dy¬ 
namics in the nucleus m , etc. Here the visco-elastic re¬ 
sponse of the system is often invoked as a physical mech¬ 
anism at hand to generate sub-diffusion (0 < a < 1), 
where the noise r]{s) and the memory kernel fj,(t) are re¬ 
lated through the fluctuation-dissipation theorem (FDT) 

fi{t - s)kBT = {r]{t)r]{s)) (3) 

with ksT being the thermal energy (see, for instance (l^ . 
fl^l and references therein). In this case, the equation 
of motion describing the process is called the fractional 
Langevin equation, which can be written in the integral 
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form as 

= y ds ti{t- s)f{s) -\-r]{t), (4) 

that is the generalized Langevin equation (GLE) with a 
power-law memory kernel, where f{t) is an external force 
acting on the system. 

In this paper, we drive such a sub-diffusive walker by 
an external force and consider its fluctuating dynam¬ 
ics. Such a situation would be related to the active 
process in the biopolymer, such as polymer transloca¬ 
tion driven by voltage drop JMl, rotational dynam¬ 
ics of entangled polymers M. and chromosome segre¬ 
gation during cell division [20l 1^ . The quantities of 
interest are the time evolutions of the average displace¬ 
ment (AD) {x{t)) and the variance in displacement (VD) 
(Aa:^(t)) (= MSD—(x(t))^), where Ax(t) = x{t) — (x{t)). 
When the memory kernel is independent of the driv¬ 
ing force, an answer would be rather straightforward; 
the linear response theory suggests the anomalous drift 
{x{t)) ~ t“/ and the fluctuation Ax{t) around the aver¬ 
age drift again becomes the fBm whose VD is given by 
{Ax'^{t)) ^ with the same anomalous exponent a as 
the undriven fBm. Then, the ratio (Aa:^(t)) / {x{t)) be¬ 
comes constant in time, which provide a way to evaluate 
the driving force from the time series of the trajectory 
(see Eq. (1321) below). The natural question here is how 
such a simple and potentially useful result should be al¬ 
tered (or not) beyond the linear response domain. In¬ 
deed, the nonlinearity in the viscoelastic response would 
be pertinent to most of soft materials and biopolymers, 
where anomalous diffusion is commonplace. The mem¬ 
ory kernel /i/(t) and the noise ?7/(t), then, become force 
dependent. 

As a model system, we analyze the motion of a tagged 
segment in a long polymer chain. This is indeed one of 
the paradigms for the sub-diffusion, where the GLE with 
a power-law memory kernel can be derived from a micro¬ 
scopic polymer model 0, HI]- In our model, the nonlin¬ 
earity arises from the self-avoidance (SA) and hydrody¬ 
namic interactions (His) between different segments in 
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FIG. 1: (Color Online) Schematic representation of driven 
dynamics in polymer stretching in (a) near equilibrium and 
(b) strongly driven dynamics. 


the polymer. Our main purposes are (i) to determine 
the nonlinear memory kernel /i/(t) and (ii) to under¬ 
stand the nature of the fluctuation of the driven tagged 
segment determined by r]f{t). Though our main inter¬ 
est is in the nonlinear (strong force) regime, we also in¬ 
clude the discussion for the linear response (weak force) 
regime. This not only makes the argument comprehen¬ 
sive, but also counterpoints the qualitative difference be¬ 
tween these two regimes. 

In Sec. II, we construct scaling theory for the anoma¬ 
lous dynamics of the tagged segment based on the notion 
of time-dependent friction. This will be done first for the 
weak driving regime (Sec. II A), and then for the strong 
driving regime (Sec. II B), where the response becomes 
generally nonlinear. In Sec. Ill, we present molecular 
dynamics (MD) simulation results in the strong driving 
regime, which are discussed in light of the scaling the¬ 
ory. The result indicates that the analogous relation to 
Eq. © 


Hf{t-s)kBT = {T]f{t)r]f{s)) (5) 

does not hold, and the fluctuation Ax{t) is not neces¬ 
sarily described by fBm. To get better insights into the 
anomalous driven dynamics, we proceed, in Sec. IV, the 
analysis of modes in a polymer. This enables us to derive 
the memory kernel from a microscopic polymer model. 
Again, we first give a paradigmatic treatment in the weak 
force regime (Sec. IV A, B), on which we develop an ap¬ 
proximate argument suitably modified to the strong force 
regime (Sec. IV C). In Sec. V, we conclude this study. 


II. SCALING THEORY 

The polymer consists of N segments, each of which is 
characterized by its size a and friction coefficient 7 . We 
prepare the polymer in its thermal equilibrium state, so 


that its average spatial size is 

R ~ aiV'^. ( 6 ) 

The Flory exponent n characterizes the equilibrium coil 
size; 1 / = 1/2 for the ideal chain; n = 3/4 or ~ 0.588 
in two-dimensional (2D) or 3D spaces, respectively, for 
the chain with SA . In addition, we write the longest 
relaxation time r of the polymer as 



where Tq ~ ja^/ksT is the segmental time scale, and 
the so-called dynamic exponent takes z = 3 in the His 
dominated case (non-draining polymer) or 2 : = 2 -|- (l/zz) 
for the free-draining polymer [23|. 


A. Time-dependent friction 


Suppose that we switch on a constant external force / 
acting on the end segment, hereafter called tagged seg¬ 
ment, at t = 0. We set its initial position as origin 
a;(0) = 0 (see Fig. 1(a)). As time goes on, the motion 
of the segment creates the tension, which gets transmit¬ 
ted along the polymer [TJ, . This leads to the time- 
dependent friction associated with the growing section 
of the correlated motion, thus gives rise to the memory 
effect to the motion. Below, we construct the scaling 
form of the memory kernel from the force balance argu¬ 
ment during the process. But the way the above scenario 
comes into effect depends on whether the motion of the 
segment is dominated by the thermal fluctuation or the 
driving force. We are thus led to distinguish weak and 
strong force regimes. For simplicity, we look at the posi¬ 
tion x{t) of the segment in the direction along the force 
only; the motion in other directions is just an unforced 
fBm. 

Weak force regime —. The force magnitude / is weak 
enough / < ksT/R so that the polymer shape is kept 
in equilibrium coil. In this weak force (near equilibrium) 
regime, the motion of the segment is essentially relax- 
ational in thermal fluctuation. Let r*{t) be the distance, 
i.e., the root VD (or root MSD), that the segment trav¬ 
els during the time interval t. In the time window t < t, 
the polymer as a whole has no time to relax. This im¬ 
plies that the restoring force ~ kBT/r*{t) acts to the 
tagged segment, and at the same time, the tension due 
to the motion is transmitted up to n*{t) ~ (r*(t)/a)^/‘" 
segments apart along the chain. These n*(t) segments 
would take part in the motion of the tagged segment, 
thus, contribute to the friction 7 (f). Therefore, the force 
balance equation reads 


lit) 


dr*{t) 

dt 


kBT 

r*{ty 


( 8 ) 


where the time-dependent friction is given by 


lit) - 1 


yt) 


z-2 


( 9 ) 
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as inspected from Eq. 0. Solving the differential equa¬ 
tion we obtain 

W(- \J (Aa;2(t))) ~ a (10) 

, ( 11 ) 


where we assume the equilibrium formula ® is valid up 
to the length scale which adds up in the larger scale. 
The force balance equation thus reads 

the solution of which is 


where Aa;(t) = x{t) — {x{t)). We thus identify the anoma¬ 
lous diffusion exponent a = whose physical origin 
is the memory effect associated with the tension trans¬ 
mission. Superimposed on this diffusion is the (small) 
drift {x{t)), which can be obtained from the relation 
'y{t)d{x{t))/dt = / as; 

where ^{t) is determined from Eqs. ® and (fTOll . Com¬ 
paring Eqs. (HQD and m, one finds that the drift is in¬ 
deed negligible as long as t <C r/o with 


TfO — To 



(13) 






ToJ 

1/2 


To) yksT J 


fa 

ksTj 


(17) 

(18) 


In addition to the anomalous drift exponent 1/2, there 
arises the force dependence with the characteristic expo¬ 
nent {z — 2)/2, i.e., the res pon se is nonlinear except for 
the Rouse model z = A . 

Superimposed on this drift is the (small) diffusion 
the property of which is to be unveiled. One 
might estimate it from the analogous relation as Eq. 
'Yf{t)dy/{Ax‘^(t))/dt ~ ksT/y/{Ax^Jt)), which yields a 
conjecture 


Given the condition / < ksT/R, we hnd Tfo > r, ensur¬ 
ing the fluctuation dominance in the weak force regime. 

Strong force regime —. Now the force is strong enough 
(/ > ksT/R) to make the drift dominant over the diffu¬ 
sion for the tagged segment motion. At the same time, 
the polymer is stretched in the pulled direction. The re¬ 
sultant steady-state conformation can be pictured as a 
succession of blobs with growing size [2^ . Eor simplicity 
of the argument, we here adopt the mono-block approx¬ 
imation |26l |. i.e., the size ^ of the blobs are uniform and 
given by^ ~ fcsT//, which suffices for the scaling discus¬ 
sion (see Appendix A and e.g. Ref for the argument 
with spatial inhomogeneity). The stretched length of the 
polymer is estimated as 


R|| 



~ a 



N, 


(14) 


where g ~ is the number of segments inside the 

blob. 

Again we switch on the force at f = 0, before which 
the polymer assumes an equilibrium state at rest. In 
the strong force regime, the driving force overwhelms 
the thermal fluctuation, thus governs the motion of the 
tagged segment. Therefore, we set r*{t) to be a drift 
distance traveled by the segment during the time inter¬ 
val t. The polymer as a whole has no time to react to 
the force, and the tension gets transmitted only up to 
n*[t) ~ lr*{t)/a){fa/kBT)^''~^'>/‘' segments apart from 
the tagged segment (e.g. Eq. (HI). These segments fol¬ 
low the driving force, hence, the time-dependent friction 
grows as 


fy(Ax2(f)) ~ a 



fa 

k^J 


(19) 


where 7 /(t) is determined from Eqs. (fTKI) and (ITTI) . 

Remarks —. 

(i) The effect of the force is a weak perturbation in the 
scale smaller than f ~ ksT/f. The corresponding time 
scale is Tfo — To(fya)^ given in Eq. (fT^ . which signals 
the onset time of the strong force regime. In the time 
range t < Tfo, the weak force regime applies [ 13 . 

(ii) The terminal time of the driven process is given by 
the condition r*(r/) ~ i?||; 


Tf ~ TqN'^ 




( 20 ) 


At this time, the tension caused by the external force 
reaches to the other chain end. This crossovers to the 
equilibrium formula (O at / —;► kgT/R. 


B. Generalized Langevin Equation 

The above scaling arguments for the memory effect 
can be generalized to the case of the arbitrary protocol 
of the time-dependent driving force f{t). This leads to 
the GLE given in Eq. (|1|) for the motion of the tagged 
segment . The equivalent expression is 

0 = -/ ds r(t--h/(t)-ba;(t), (21) 

J-oo as 


7 /( 1 ) - 7 


aj g ^ \ a J 



where the friction kernel r(/) and the noise uj{t) are re¬ 
lated to the mobility kernel fi{t) and the noise ?/(/) in 
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Eq. dU as f(z)/^z) = 1 and w(z) = t{z)fi{z) in the 
Laplace domain Q- 

lEeafc force regime —. In our protocol, we switch on a 
constant force at t = 0 , i.e., f{t) = / 0 (t) with 0 (t) being 
the Heaviside step function. The time-dependent friction 
is then 7 (t) = [/J ds Using Eqs. (O and (ITUl) . 

we obtain 


- 

m ^ 



( 22 ) 

(23) 


where the minus sign in fi{t) comes from the fact 2 — z < 
0 in practice. This reflects the viscoelastic response of 
the tagged segment leading to the sub-diffusion. Note 
that such a power-law regime should be viewed as an 
intermediate asymptotics valid in the time range tq 
t <C T (See Sec. lIVI) . We assume that the memory kernels 
(/i and r) are related to the noises {rj and to) through 
FDT given by Eq. ([3]); the equivalent expression is r(t — 
s)kBT = (w(t)w(s)). The AD and VD can be calculated 
as 


{x{t)) = ds dufi{s-u)f 

(Aa;^(t)) = J ds J du {r]{s)ri{u)) ^^ . (25) 


in agreement with Eqs. dni) and OT, respectively. 

Strong force regime —. Assuming the same line of argu¬ 
ment as in the weak force regime, we obtain the estimate 
of the nonlinear memory kernel 


r/(i) 


1 

7^0 


-3/2 


T 

To 


To 
- 1/2 


fa 

ksT 


(z/2)-2 


fa 

ksT 


2-(z/2) 


(26) 

(27) 


One can easily check that this yields the drift scaling, 
which is in accord with Eq. dni- The fluctuation around 
this drift is subtle, however. In Ref. [1^, it was assumed 
that the property of the noise rjj (or w/) is encoded in the 
kernel through the relation ([S]). With this naive assump¬ 
tion, the solution of the GLE leads to the VD scaling, 
which is in accord with Eq. dn}. We repeat once more 
that this estimate needs to be checked. 

Summary —. The result and conjecture obtained so 
far based on the scaling argument are summarized as 
follows: The leading component in the dynamics of the 
tagged segment is fluctuation or drift in weak or strong 
force regime, respectively. In either case, the motion of 
the tagged segment creates the tension, which gets trans¬ 
mitted along the chain with characteristic dynamics in 
respective regime. The resultant anomalous dynamics 
can be expressed using the memory kernel as 


(Ax^it)) 


keT 

r(t) 


[weak force regime] (28) 


{x{t)) 


f 

r/(i) 


[strong force regime] (29) 


The FDT ([3]) or its analogous relation ([5]) then implies 


{x{t)) 

{Ax^{t)) 


J_ 

m 

keT 

r/(i) 


[weak force regime] (30) 
[strong force regime] (31) 


for the minor component in the motion. In addition, it 
is worth pointing out that GLE formalism with FDT ([Sj) 
or its analogue dS]) predicts the quantitative relation be¬ 
tween the drift and the fluctuation 


(Ax^(t)) = ^(x(t)) (32) 


regardless of the specific form of memory kernel. Very re¬ 
cently, this relation has been used to evaluate the driving 
force of the bacterial chromosome segregation in vivo . 


III. MOLECULAR DYNAMICS SIMULATIONS 


In this section, we perform MD simulation to verify the 
scaling predictions in Sec. II, in particular those Eqs. dH 
and (1311) in strong force regime. In simulations, equation 
of motion for each segment is 

where /(I ^ N) is bead indices, m and 7 are mass and 
frictional coefficient of a bead, Ci{t) is a Gaussian white 
noise with mean zero and the variance = 

2^kBTdijS{t — t')l with I being a unit matrix, and 
is the unit vector directed to x-axis. His are ignored 
for simplicity (free draining). The total potential U = 
UpE + Urep consists of the finitely extensible nonlinear 
elastic potential; 


UpE = 


UpE 


N-l 

( 2 a )2 log 

2=1 


b+i^\(34) 

(2a)2 ^ 


and the repulsive potential between different beads; 

^12 

i<j ' * 

= 0 for \xi — Xj\ > R. 

We set e = kBT, Gpe = lOfe^T/a^, 7 = (jnkBTj 
and the unit time is tmd = mj^ = {ma^/kBT)^^"^ = 
'JO?/kBT. The total number of beads is A = 100, for 
which the equilibrium size is i? ~ 6 . 8 a in 3D and ~ 9.5a 
in 2D. Initial conditions are picked up from equilibrium 
conhgurations, and the force applied at A-th bead is 
switched on at t = 0. We choose the time step for inte¬ 
grating equation of motion as 5t = 0.005rMD- 
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FIG. 2: (Color Online) Dynamics of the tagged segment 
from MD simulations. Left or right column shows results in 
3D or 2D, respectively, (a) Time evolution of average drift 
and VD at fa/ksT = 0.3. Force dependence of (b) aver¬ 
age drift and (c) VD at time I/tmo = 5 x 10^. All plots 
are displayed in a double-logarithmic scale. Inset triangle 
slopes indicate the theoretical scaling exponents speculated 
by Eqs. (|27ll.(|29l).([3T1l. 


Figure 2 (a) shows the time evolution of AD {x{t)) 
and the VD {Ax^{t)) of the pulled bead by the force / = 
OAksT/a in 3D (left) and 2D (right). Both quantities 
exhibit a slope close to 0.5 in agreement with scaling 
predictions (II3 and m- After long time t/rMO > 5 x 
10^ (10^) in 3D (2D), the dynamics becomes normal, i.e., 
{x{t)) ^ {Ax‘^{t)) ~ t, which corresponds to the center 
of mass mode. This time is interpreted as the tension 
propagation time r/ Eq. (1^ . at which the effect of the 
driving force reaches to the other side of the chain end, 
and the chain settles in the steady-state as a whole. 

Figure 2 (b) shows the force dependence of the drift 
{x{t)). The data is taken at time I/tmo = 5 x 10^ in the 
anomalous dynamics regime (see Fig. 2 (a)). Recalling 
the dynamical exponent z = 2 -\- v~^ in the free draining 
case, the results are in good agreement with the predicted 
force exponent from Eq. (flTl) . that is {2v)~^ ~ 0.85 (0.67) 
in 3D (2D). 

Figure 2 (c) shows the force dependence of the VD 




th 


MD 


t/z 


MD 


FIG. 3: (Color Online) VD-AD ratio 

f (^Ax^{t)'j /{ksT {x{t))) obtained from MD simulations 
in 3D (left) and 2D (right) as a function of time. 


{Ax‘^{t)). Again, the data is taken at time t/xMO = 
5 X 10^. The predicted force exponent from Eq. m is 
(2z/)-i - 1 ~ -0.15 (-0.33) in 3D (2D). We find non- 
negligible deviations from the conjecture m, which is 
more evident in 2D. 

Next, Fig. 3 shows the normalized VD-AD ratio 
/(Aa;^(t ))/{ksT {x{t))). Recall that the GLE with the 
relation m predicts this ratio to be time independent 
and exactly 2 (eq. (15^ 1. While in 3D, this relation is 
satisfied reasonably well over a wide time range, appre¬ 
ciable deviations are found in 2D, which gets larger as 
the force becomes larger. 

The above results on the fluctuation Ax{t) of the 
tagged segment around the average drift suggests that, 
against a naive expectation, it cannot be described as 
the fBm in the strong force regime. What causes it? To 
answer this, we shall take a microscopic polymer model, 
and attempt to derive the memory kernel and GLE based 
on the analysis of modes in the polymer chain. 


IV. MODE ANALYSIS 

We first provide exact calculation for the Rouse model 
and deliver some remarks. We then approximately incor¬ 
porate the nonlinearity due to SA and His into the mode 
equation by making the effective friction and spring con¬ 
stants mode-number dependent (Sec. lIVBl) . Based on the 
solid framework on the weak force regime, we develop 
an approximate treatment in the strong force regime in 

Sec.ICT 


A. Rouse model 

The polymer is modeled by A -|- 1 connected beads. 
These beads have no excluded volume, linked by har¬ 
monic springs in series (the root-mean-square length of 
each spring is a), and move in a viscous fluid by being 
kicked by thermal noise and external force. As in Sec.llll 
we again focus on the direction of the force only, and sup¬ 
press the vector notations. In the limit where the bead 
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labeling index n £ [0,1V] is made a continuous variable, 
the equation of motion reads 

where the friction coefficient for a bead 7 and the spring 
constant k = iksT /a^ (in 3D) defines a microscopic time 
scale To = 7 /A:, and fn is the external force acting on n- 
th bead. Open boundary conditions are imposed at both 
chain ends for linear polymers 

dxji dxji 

dn n=o dn 

The random forces Cn(t) acting independently on indi¬ 
vidual beads are Gaussian white noise with zero mean, 
whose correlation obeys the FDT; 

(Cn(A)Cm(s)) = 2jkBTS{n - m)6{t - s). (37) 

In the Rouse model, the noise and the external force 
acting on some segment affect the motion of other seg¬ 
ments through the elastic connectivity. We may thus 
expect that the above FDT (1571) at the individual seg¬ 
ment level could be coordinated to generate a relation 
between the fluctuation and the response at the level of 
collective dynamics of the entire chain. Below, we shall 
derive such a relation 

{Vnit)Vm{s)) = kBTXnm{t,s) (38) 

with the concrete functional form of Xnm(t, s) and r]n(t). 
Here the response function Xnm(t,s) = 5{xn(t))/5fm(,s) 
describes the change in the average velocity of n-th seg¬ 
ment at time t caused by the force that acted on m-th 
segment at time s(< t), and r/n{t) is a correlated Gaus¬ 
sian noise with zero mean {rjn{t)) = 0 acting on n-th 
segment. The FDT (1551) thus indicates that the cross 
correlation of the noise has a long time memory, which 
is related to the collective response of the segment. It 
includes Eq. ([51) as a special case of the self-response and 
correlation n = m. 


To this end, we analyze the Rouse equation (1551) in 
terms of the normal coordinate Xp(t) 

pN 

— / dn hp fi^nid^ 

Jo 

(39) 

with 


1 /TrpnN 

Kn- jyCOS^ ^ j. 

(40) 

Its inverse transform is 


Xn{i) = ^ ^ 

p >0 

(41) 

where 


2 <=«(T) 

1 + dpo ■ 

(42) 


The normal modes obey the following equation of the 
overdamped harmonic oscillator type 

= -kpXp{t) + Zp{t) + Fp{t), (43) 

where the spring and the friction constants kp^ jp de¬ 
fine the relaxation rate of the p-th mode kp/^p = 
{k/'y){TTp/N)'^, and Zp = J dn {'yp/j)hp^nCn{t) is the 
noise in the mode space. The external force can be 
arbitrary, but for our present purpose, we manipulate 
a particular segment (labeled by the index m), i.e., 
fnit) = fmit)5(n — m), which is distributed in the 
mode space according to Fp(t) = J dn {'jp/'y)hp^nfn(t) = 
{lpll)hp,mfm(t). There is no restoring force for p = 0 
mode, which corresponds to the motion of the center of 
mass. 

It is useful to set as Xp = 27V7/(1 -|- ^po), thus kp = 
2A:(7rp)^/iV, so that the FDT in the mode space takes a 
familiar form: {Zp(t)Zg{s)) = 2^pkBTSpq6(t — s). Equa¬ 
tion (|55|) is solved as 

Xpit) = — [ ds e-('='>/>h‘-^)(Zp(s) + Ep(s)) 

7p Jto 

(44) 

where Xp(to) is the initial condition for the mode p. We 
assume to — oo so that the system is in a stationary 
state before we apply the external force. By direct calcu¬ 
lation, one can check the following FDT (see Appendix 

B) 

Cp{t, s) = kBTxpit, s), (45) 

where the velocity correlation function Cp{t,s) 
and the response function x(t, s) are defined as 
{AXp{t)AXq{s)) = 6pqCp{t,s) and Xp(t, s) = 
S{Xp{t))/SFp{s), respectively. 

Upon time derivative of Eq. (I44p and transforming it 
into the real coordinate using Eq. dm , one can express 
the time derivative of the position of n-th segment in the 
following form; 

—^ = J ds Xnfnit,s)f^{s)+ r]n(t), (46) 

where 

X-nm '^) E Xp{t,s)hlnhlrn ( 47 ) 

p>0 

and 

Vn{t) = E / Xp{t,s)Zp{s) (48) 

p^od-oo 

are, respectively, interpreted as the velocity response 
function and the noise (The functional form of Xp(t, s) 
is given in Eq. (1551) 1. The latter is identified as 
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the velocity fluctuation r 7 „(t) = Ain(t) = Xnit) — 
{xnit)) of n-th segment. The fluctuation-response re¬ 
lation (IMl) can be most easily verified by decompos¬ 
ing the response and the correlation functions into 
modes: Eq. (gT]) and (??n(t)??m(s)) = {Axn{t)Axmis)) = 
J2p>o Cp{t, The FDT (|45|) in the mode space 

then leads to Eq. (1551) . 

Memory kernel —. Suppose we tag a particular seg¬ 
ment (labeled by m). The force /(t) is applied only 
to that m-th segment, and we track its stochastic mo¬ 
tion. The information of interest is contained in the self¬ 
response s). In this context of the single segment 

tracking analysis, positions of other segments except for 
m are inaccessible to observations. One can therefore 
omit the label index and write the equation of motion of 
the tagged segment as 

= J ds s)f{s) + r]{t) (49) 

where, as verified above, the FDT ([3]) holds. The mobil¬ 
ity kernel Xmm(t, 0)) = + no{t) + finit) 

composes of three terms according to Eq. (H7)) : 


the Rouse polymer. The minus sign here is a hallmark 
of the viscoelastic response inherent in the system with 
elastic connectivity (see, for instance, a similar analysis 
for polymerized membrane [1^). 

(hi) In Eq. (l52l) . the mode 


P*{t) 


V/ 


(53) 


has the largest contribution at time t. The corresponding 
number of segments is n*{t) ~ N/p*{t) ~ (t/ro)^^^. This 
agrees with our scaling estimate for the tension front in 
the weak force regime (Eq. (|lip with z = A and i/ = 1/2). 
The physics behind this agreement is the following; The 
effect from the larger scale beyond the tension front n* (t) 
is irrelevant, or only marginal. Therefore, even if we 
neglect it, i.e., by shifting the lower bound pi,b in the 
summation pri, = 1 to pi,b = P*(t), one should get a 
qualitatively correct result with the proper exponent, 

(iv) The following sum rule 

dt [poit) -I- PAiit)] = 0 (54) 


,r\t) 

Po{t) 

PM{t) 


-m = ^d(t) 

70 -/V 7 


N 




N 




kj 

p=l ’P 

1 1 
40r7ro V'^o 


- 3/2 


(50) 

(51) 

(52) 


(to < t < r), 


where r = 71 /fci ~ tqN'^ is the terminal (Rouse) time. In 
the above, the upper bound in the summation over p is set 
reflecting the original discrete nature of the model with N 
degrees of freedom, and we replace cos^ {pirn/N) by the 
average 1/2. The last near-equality in Eq. (1551) is valid in 
the time window tq <IC t <C t, where the summation over 
p can be replaced by the Gaussian integral. For longer 
time t > T, only the p = 1 mode is relevant, thus, the 
memory decays exponentially. 

Remarks —. 

(i) In the very short time scale {t ~ tq), the seg¬ 
ment is unaware of the connectivity, and exhibits a usual 
viscous response with the segment friction coefficient 7 
(Eq. (I5T]) 1. In the time scale coarser than r, the response 
again becomes viscous, but now with the much larger 
friction coefficient qfV, i.e., the center of mass mode 
(Eq. (1551) 1: see the sum rule Eq. (IMl) below. 

(ii) In the intermediate time scale (tq t <S; r), the 
last term pM{t) (Eq. (1551) 1 dominates the dynamics of 
the tagged segment. Noting z = 4 for a Rouse model, 
the result dsn agrees with the scaling analysis Eq. (1551) 
in Sec. (HD. This term represents a memory effect, which 
arises from the superposition of the internal modes in 


may hold for any physical system, which indicates that 
in the long time limit t ^ r, all the internal modes relax, 
and only the center of mass mode remains. 

(v) In certain visco-elastic environments, the segment 

response itself could contain the memory effect. Then, 
one may think of the visco-elastic Rouse model, where the 
viscous friction term jXn in the Rouse equation (1551) i s re¬ 
placed with the integral kernel f ds ri(t—s)xn(s) isn. 
In the case of the power-low memory kernel ri(t) ^ 
with 0 < «! < 1 , the exponential relaxation in vis¬ 
cous Rouse model (Eq. (04])) is generalized to the non¬ 
exponential one described by the generalized Mittag- 
Leffler function. This results in the memory kernel in the 
tagged segment dynamics pM{t) ~ hence 

the anomalous exponent a = ai /2 for the tagged seg¬ 
ment diffusion. Such a visco-elastic Rouse model has 
been proposed to analyze the sub-diffusion of bacterial 
chromosomal loci . The usual viscous result cor¬ 

responds to the limit ai —^ 1. The relation between expo¬ 
nents for the single segment exponent ai and the tagged 
segment one a with a factor 2 is a general consequence 
of the Rouse model. 

(vi) The Rouse model is valid as long as / < ksT/a. 
For larger force, the chain section close to the pulled site 
becomes highly stretched, forming a “stem” . For the 
prescription and the scaling analysis in such a situation, 
see Ref. M- 


B. Self-avoidance and hydrodynamic interactions 

In many of practical situations, one of or both of 
these effects become important. These interactions are 
non-local, and conformation dependent, hence, make the 
equations of motion highly nonlinear, which prevent the 






rigorous analysis based on the mode expansion. Nev¬ 
ertheless, one can gain physically appearing insights in 
terms of approximate mode analysis. The pre-averaging 
approximation provides a way to treat the His in terms 
of modes, in which the conformation-dependent mobility 
tensor is averaged using the equilibrium segment distri¬ 
bution . This yields for the effective friction constant 
for the mode p 


^ (p = 0) 

\ 7p(Ar/p)^(^-2) (p ^ 0) ^ ’ 

Note that the free-draining polymer z = 2 + {l/i') bears 
no relation to the pre-averaging; we then recover 7 p for 
the Rouse model. 

In a similar level, the self-avoidance (SA) can be 
treated by employing the linearization (Gaussian) ap¬ 
proximation, which alters the spring constant for the 
mode p as 

kp ~ kp{p/N)^''. (56) 

The validity of this form as well as a high degree of sta¬ 
tistical independence of different modes has been numer¬ 
ically demonstrated in Ref. [s^. Note that for the ideal 
polymer v = 1/2, we recover kp for the Rouse model. 

Let us analyze the mode equation (1431) with Eqs. dSS]) 
and (1561) . Note that the terminal time is now given by 
T = 7 i/A:i ~ (see Eq. (O). The mobility kernel is 

again decomposed as pi{t) = + p^it) + pnit)- 

While the segment instant response po{t) is essentially 
unchanged from the Rouse model (Eq. (I51|) l. the N de¬ 
pendence of the center of mass response is modified as 
~ In addition, the memory 

kernel is evaluated as 


p—1 ' 


1 / t 


ITo \To 


-{2-2z~ 


(to < t < r) (57) 


The last near-equality is valid in the intermediate time 
scale, where the summation is evaluated as the integral 
using the formula^ /o°°= T[b/6)/{9a^^^) 
for a,b,9 > 0. The result agrees with our scaling argu¬ 
ment (see Eq. (12^ 1. and the tagged segment dynamics 
in this time scale is a fBm with the anomalous exponent 
a = 2lz. 

Remarks —. 

(i) At time t, the mode p*{t) ~ has the 

largest contribution. The corresponding number of seg¬ 
ments n*{t) ~ N/p*(t) agrees with our scaling argument 
Eq. dill) for the dynamics of tension front. 


^ Tlie symbol r(-) here is used as the gamma function r( 2 ) = 
u^~^e~'^du in this integral formula. 


(ii) For the present description to be valid, at least 
qualitatively, the condition / < ksT/{aN'^) is required. 
For stronger force, the conformation of the polymer 
is markedly deviated from the equilibrium distribution, 
which invalidates the assumption used to evaluate effec¬ 
tive friction and spring constants in Eqs. (ESI) and (EH). 
This is contrasted to the Rouse model case, for which 
the condition is much weaker, associated to the bond 
stretching (See remark (vi) in Sec. IIV Al) . but not the 
conformation. 

In Sec. IIV Cl we aim at constructing an effective de¬ 
scription, which allows us to analyze the fluctuating 
driven dynamics in larger force regime / ^ ksT/{aN'^) 
even with SA and His. 


C. Driven dynamics 

Suppose we start applying a constant strong force 
/> ksT/{aN'^) to the chain end (V-th segment) at time 
t = 0 (see Fig. 1(b)). Before that moment (t < 0), the 
polymer assumes an equilibrium conformation. We are 
interested in the motion of that pulled segment. The dy¬ 
namics is nonlinear with SA or His and nonequilibrium 
in strong force regime. To analyze the average dynamics, 
the following nonlinear diffusion-type equation (called p- 
Laplacian diffusion equation) has been proposed [S^ [Sg : 



with p = (z — 2)z//(l — v) and the segment diffusion 
coefficient Dq ~ ksT/^. This equation is derived based 
on the force balance argument for the chain of blobs, 
and can be thought of as a nonlinear extension of Rouse 
model (see Appendix A for the derivation). 

Equation (1551) would be a useful starting point to an¬ 
alyze the stochastic dynamics of the polymer stretching 
in terms of the collective mode in the chain. Since, by 
construction, this is expected to provide a reasonable de¬ 
scription on the average dynamics, one may add a ran¬ 
dom force Cn of zero mean to get a nonlinear Langevin 
equation 


= kn^^ + Cnit) + fn, (59) 

where and 7 „ are given by Eqs. (l88)) and (l89l) in Ap¬ 
pendix A. Within the mono-block approximation, one 
can linearize it 

with the force-dependent spring and friction coefficients 
Eqs. EQl), ED)- The random forces Cn independently act¬ 
ing on individual segments are assumed to be Gaussian 

/ j:\ 

white noise with zero mean (Cn {t)) = 0, whose corre¬ 
lation obeys the FDT; {Cn\t)(m\s)) = 2j^^'>kBTS{n — 
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m) 5{t — s). Equation (l60l) reduces to the Rouse model 
(eq. ( 1551 ) ') when v = 1/2 and z = 2 + {1/v) = 4. Oth¬ 
erwise, the SA or His result in the nonlinear response. 
With the force free boundary condition (Eq. ( 1551) 1 and ex¬ 
plicit inclusion of the external force fn{t) = fQ{t)S{N — 

n) acting on the end segment n = N, one can follow the 
analysis developed in Secs. IIV Al - IIVBI 

In analyzing the dynamics of nonlinear response, one 
has to be aware of the change in the mode spectrum, 
i.e., (fcp, 7 p) due to the external force. 

This effect may be treated in the following way. After 
the force is switched on at t = 0, the equilibrium mode 
would persist during the induction time (see remark (i) in 
Sec. Ill Ap . Equation (|43)l with Eqs. (1551) and (|56ll would 
be thus valid up to t = r/o- At t > t/o, the effect of the 
driving dominates the mode dynamics, and the spring 
and the friction constants become altered to those in the 
strong force regime. Thus, making these constants time 
dependent, the equation of motion in mode space be¬ 
comes 

7; = -K^t)Xp{t) + z;{t) + Fp(t), (61) 


where the equilibrium mode structure 7 p(t) = 7 p, k*{t) = 
kp, Z*{t) = Zp{t) persists only up to t < r/o (Eqs. (l55l) 


and (l56]) l. At t > t/q, these switch to the stretched mode 
7p(0 = lp\ k*{t) = kl^\ Z*{t) = Z^^\t) with 


^U) ^ f {p < Pf) 

P 1 ip{N /(p > Pf) 


(62) 


r(/) 

h,p 


Yf)p^/N {p < Pf) 
kpijp/N)'^'' {p>pf) 


(63) 


and the Gaussian white noise Zf (t) with zero mean and 
the correlation {Zp^\t) z\^\s)) = 2^l^^kBTSpq6{t — s). 
The tension propagation time given in Eq. dSni) can be 
identified as the slowest relaxation time in the stretched 
mode Tf = /k[^\ Here, the characteristic mode num¬ 

ber is introduced 


Pf = 2N 



N 

9 


(64) 


in such a way that the effect of the force is negligible for 
modes with p > Pf{<^ a{2N/p)'' < ksT/f), hence, the 
friction and spring constants are given by those in weak 
force regime (Eqs. (ISSl) and (l56l) . respectively). Notice 

(p) 

that our construction assures continuity for both 7 / and 
k^/^ across p/. 

The solution of Eq. dUD for t > T/o takes the same 
form as Eq. (0 with the replacement (fcp, 7 p,to) —>■ 
{kp \jp ,Tfo), where the “initial” condition is given by 

1 /•’■/o fcp , , 

XpiTfo) = — / ds e Zp{s) 

7 p J—00 

+ ^ . (65) 



102 102 


lO'* 


10= 



lO'* 10= 

th^ 


FIG. 4: (Color Online) VD-AD ratio f (Ax^) /{kgT {x)) 
obtained from theory (Eq. (1731) 1 in 3D (left) and 2D (right) 
as a function of time. The calculation was carried out with 
N = 100. Note that there are uncertainties in the precise val¬ 
ues in the peak heights due to the scaling estimate of spring 
constants. In the above plots, we set all the numerical con¬ 
stants (including (/i^ jy)^) unity. 


Crucially, this noise from the “initial” condition adds an 
extra contribution to the noise 77 /(t) in the motion of the 
tagged segment, which leads to deviation from the rela¬ 
tion (IS|) due to the SA effect as will be discussed below. 

The motion of the tagged segment is described by 
Eq. (|55)) in the weak force regime at t < t/q. To take 
account of the changes in spring and friction constants 
at t > T/o, we need to modify it as 


dx{t) 

dt 


/ ds pf{t,s;Tfo)fis)+r]f{t;Tfo) (66) 

J —00 


with the mobility kernel 

Pf{t,s;Tfo) = ^Xp,/(t,s;T/o)(/i^_^)^ 

p>0 

and the noise 


(67) 


VfitlTfo) = I dsXpj{t,s;Tfo)Z*{s) 

p>0 




k;{t) 




kp{s) 




0 %[s) 


where Xpj(t, s;t/o) is the velocity response function for 
the mode p calculated in Appendix B. 

Memory kernel —. The mobility kernel /i/(t,s; t/q) in 
the intermediate time scale Tfo <C s < t <C t/ is sta¬ 
tionary, and dominated by the memory effect due to the 
connectivity, that is pf{t, s;Tfo) ~ PMj{t — s) with 


Ml 


■vp 


p>l WP ) 

Pf 


= -E: 


1 


p 


^ -fNTfo \pf J yksT 




e p 

2 / r \ Z — 2—V' 


(^p,Af)^ 


^ jNto \nJ ^ pX’ 


p=Pf 


(69) 
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In the time window r/o <S^ t ^ Tf, the second term is 
negligible and the summation in the first term can be 
approximated by the Gaussian integral. This calculation 
leads to our scaling estimate (1^ in Sec. HIl For longer 
time t > Tf, the memory decays exponentially, and only 
the center of mass (p = 0) mode remains. It is character¬ 
ized by the viscous response with the friction coefficient 

Fluctuation and response relation —. The correlation 
function of the velocity fluctuation is 

{vf{t',Tfo)vf{s;Tfo)) = J2c,j{t,s;rfo)ihl^r, (70) 

p>0 


where Cpj{t,s) is the correlation function in the mode 
space calculated in Appendix B. Comparing this with 
Eq. (1571) . we find 


{Vf{t',Vo)Vf{s-,Tfo)) - kBTpf{t,s;Tfo) 

f 0 {s <t < r/o) 

= < 0 (s < r/o < t\ 7 l) 

[ Ep>i C^pj{t,s]Tfo){hlj^f {Tfo <s<t) 

Here, a factor from the “initial” condition can be evalu¬ 
ated using Eq. dSIl). This leads to “zero” in the first and 
second lines, and the expression of the excess term in the 
third line in terms of the change in the spring constant 


C^‘;f(t,s;Tfo) 


f keT kBT\ 

) \ h ) 




(72) 


This term breaks the time translational invariance, and 
decays exponentially with a characteristic rate kp / 7 p . 
From this, we obtain^ 


{^x{tf) - 



2kBT 


{x{t)) 



(73) 



The deviation (right hand side in eq. (l73ll l is positive 
due to the force induced hardening kp ^ > kp (for SA 
chain). Comparing the 2D and 3D cases, the larger 
deviation is expected for 2D as fc), 2 d > ad¬ 

dition, Eq. dZSl) tells that the deviation grows as time 
evolves until r/. This means that the VD-AD ratio 
f{Ax‘^{t))/{kBT{x{t))) peaks around r/. At f > r/, 
p = 0 mode in the denominator overwhelms the inter¬ 
nal modes p > 1, leading to the recovery of the relation 


2 Begining with the solution Xp{t) of Eq. ED makes the calcula¬ 
tion easier to check Eq. (El. 


/ (x) = 2kBT. These trends are clearly seen in 
Fig. 4, where we plot the VD-AD ratio obtained from 
the above theory. Our present treatment is rather crude 
in the sense that the switching to the strong force regime 
around t ~ Tfo is treated through a discrete jump in the 
effective parameters. In reality, it would take place more 
smoothly. Nevertheless, our theoretical prediction well 
captures the essential trend in the MD simulation results 
in Fig. 3. 

Remarks —. 

(i) At time t, the mode p*it) ~ has the 

largest contribution in Eq. (1551) . This corresponds to 
the number of segments n*{t) ~ N/p*{t), which agrees 
with our scaling estimate (1751) for the tension front in the 
strong force regime. 

(ii) The mode with p < p*(t) may be unphysical for 
the stretching process as such a large scale part is not 
stretched by the force yet. However, we expect that these 
fictive modes do not alter the qualitative conclusion on 
the dynamics of tension propagation, just as the case in 
the weak force regime (see remark (iii) in Sec. IIV All . 

(iii) A rough estimate in the peak height in VD-AD ra¬ 
tio can be obtained by evaluating p = 1 mode in Eq. dZSl); 


f{Ax^{Tf)) 
kBT {x{Tf)) 


/T-//70 


2-b 


V kBT 


if) 

I 


ju-i _ r^(/)]-i 

/ ■. ' (1 - 

-1 CO (74) 


With a factor cq ~ (1 — e ^)^(/ip 7 v)^ ~ we obtain 
VD-AD ratio ~ 5.6 and ~ 2.4 in 2D and 3D cases, re¬ 
spectively, with N = 100 and fa/kBT = 1. Despite 
indefiniteness of these numerical values (see the caption 
of Fig. 4), this estimate would be useful to see the qual¬ 
itative dependence on the force and the chain length. 

(iv) The tension propagation time r/ fluctuates in each 
realization of the stretching processes, and Eq. (20) is re¬ 
garded as the average (r/) = /k[^\ In the strong 

force regime, the dominant source of the stochasticity 
comes not from the noise Zp(t) but from the initial con¬ 
formation of the polymer along which the tension prop¬ 
agates. In terms of the mode analysis, the fluctuation in 
Tf can be evaluated in the following way. 

For clarity of the argument, suppose that the force is 
strong enough / ~ kBT/a so that the induction time 
is very short r/o ~ tq. Neglecting the noise effect, 
the displacement of the center of mass {p = 0) mode 
and the slowest relaxational {p = I) mode are, respec¬ 
tively, Xo{t) - Ao(0) = Fot/j^^'^ and Xi{t) - Ai(0) = 
(Fi/fcJ'^^ —Ai(0))(I —6“^!^'*/'’'!^'). Comparing these, one 
finds that at f = (r/), the displacement in p = 0 mode 
reaches the final displacement in p = I mode on average, 
i.e., (r/) ~ 7 Q'^^Ei/(FofcJ'^^) ~ 7^'^V^i'^^- Taking account 
of the fluctuation in the latter due to the initial distribu¬ 
tion, one obtains the variance in the propagation time as 
((Ar/)^) ~ (7Q'^^/Eb)^(Ai(0)^). Evaluating the variance 















11 


ill the initial distribution using Eq. dSIl), this leads to 

i - / fa 

v«^ == j (™) 

The same result has been obtained in ref. [s^ using the 
scaling argument. 

(v) Events on the short length and time scales are de¬ 
scribed by the weak force regime. Such a range increases 
with the decrease in the force; the characteristic mode 
number and the time change from pf N and Tfo ~ tq 
at / ~ ksT/a to p/ ~ 1 and Tfo ~ r at / ~ ksT/R. 

(vi) Our theory indicates that it is the nonlineality in 
the elastic response (force-dependent spring constant), 
but not the frictional response, that is responsible for the 
non-trivial VD-AD ratio in the stretching process. The 
sign of Eq. (1751) depends on whether the system exhibits 
the stiffening or the softening under the force. 


response in the driven process, which is in a rather good 
agreement with results obtained from MD simulations. 

A recent study has made use of the VD-AD ratio of the 
labeled locus of bacterial chromosome during the segre¬ 
gation process to estimate its driving force [^. We feel 
that our present study could be a useful guide for such 
an analysis. 
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V. CONCLUDING REMARKS 


Appendix A 


It has been long known that a tagged segment in a 
polymer exhibits a sub-diffusion in the intermediate time 
scale, and its consequence ranges from the dynamical 
function of biopolymers to the rheology of polymer so¬ 
lutions. In this paper, we formulated the problem in 
terms of the mobility problem, i.e., the dynamical re¬ 
sponse of the segment after the application of external 
force, and introduced the weak and strong force regimes 
for the anomalous dynamics. 

In the weak force (equilibrium) regime, the motion in 
the intermediate time scale is dominated by the mem¬ 
ory effect, leading to a conventional fBm. We performed 
lucid and exact analysis for a Rouse model, which leads 
to a microscopic basis for the fBm. The deduced mem¬ 
ory kernel fully agrees with a simple scaling argument 
based on the physical picture of tension transmission. To¬ 
gether with the approximation scheme to include the SA 
and His, we believe that the present approach provides 
a comprehensive picture on the anomalous dynamics of 
the tagged segment in the weak force regime. 

In the strong force (driven) regime, the motion in the 
intermediate time scale is again dominated by the mem¬ 
ory effect arising from the tension transmission, but now 
the tension dynamics accompanies a large conformational 
distortion, and qualitatively different from that in the 
weak force regime. We discussed that the memory kernel 
generally becomes force dependent, from which one can 
derive the nonlinear dynamical scaling for the anomalous 
drift. Unlike the weak force regime, the fluctuation and 
the response do not necessarily satisfy a simple propor¬ 
tionality relation due to the noise from the “initial con¬ 
dition” at t = Tyo, after which the dynamics enters the 
strong force regime. This extra noise is non-stationary, 
making the fluctuating dynamics to deviate from the 
fBm. On the basis of the approximate mode analysis, 
we proposed a formula to relate the fluctuation and the 


We review the previous works on the deterministic (av¬ 
erage) dynamics. 

Coarse-grained description: a chain of blobs —. 

To discuss the dynamics in the scale larger than the 
blob size, we envision the stretched polymer as a chain of 
blobs. The blobs are labeled with the index h = 0,1, 2, • • • 
from the free end at the rear. The h-th blob comprised 
of pn segments has the spatial size £,n — By taking 
the X axis as the pulling direction, the position in of the 
center of h-th blob is 


n—1 

Xfi = Xn—Q + ^ ^ 

(76) 

h—0 


The dynamics of the chain of blobs can 

be analyzed by 

noting that the spring and the friction constants kn, ja 

of the h-th blob are given by 


r keT 

(77) 

I/’ ~ 

- .2 

7., == 7 ( 1 ) . 

(78) 


which lead to the force balance equation 


with the elastic restoring force 


fn — kn(^Xn-\-l Xfi) kn—l(Xn Xfi—l) 


A 

dh 


din 
‘ dn 


(continuum limit) (80) 


and the viscous frictional force 

r{vis) _ ~ din 

In - 


( 81 ) 
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Equation ((76l) indicates the expression for the local de¬ 
formation 


dxn _ . 

- SS 

on 

One can then write Eq. (Eni) as 


(82) 


This line of argument was used to discuss the normal 
modes of the tethered chain stretched by flow [2^. 

Mapping to the p-Laplacian equation — 

One may extrapolate the above force estimation at the 
blob scale to the segment scale in such a way that the 
elastic and viscous frictional forces acting on the n-th seg¬ 
ment are given by and 

The label index of blobs and that of segments are related 
as n = gfi'dh'. We write respective forces as 


/r> = (84) 

4'“’ = -7.^ (85) 


The spring and the friction constants 7 „ in this hne- 
grained frame can be estimated in the following way. The 
relation dn = gndh between internal coordinate before 
and after the fine-graining indicates the transformation 
rule of the local chain deformation 


dXn _ dxn 

a~ — 9n~^ 

on on 


( 86 ) 


This, together with Eq. (1821) . implies 


a 




(87) 


These considerations lead to ^ 


~ fc ( — 

a 


{l-2v)/v 


7„ ~ 7 — 

a 


z-2-{l/v) 


~ k 
dx. 


dx„ 




~ 7 


d{na) 


( 88 ) 


d{na) 


(89) 


The force balance relation = 0 can be cast 

into a so-called p-Laplacian diffusion equation given in 

Eq. dig) [Mill: 

Again, useful insights can be deduced from the mono¬ 
block approximation [^ . where the blob sizes are as¬ 
sumed to be uniform with — ksT/f (see Sec. Hill . 


® The second derivative relation also follows as d^Xnldfi^ = 
Cg? d^Xnldr? with a negative coefficient C = {yjiy — 1)) < 0. 


In this approximation, the spring and the friction con¬ 
stants Eqs. (1551) . (1551) do depend on / but not on n: 


u — h(f) ^ t. 
^ ^ 

7n = — 1 


r y2.-i)/. 
fa 

k^J 


(90) 

(91) 


Therefore, Eq. (|58)) becomes a simple linear diffusion 
equation 


/ /g Y d'^Xn 

dt ° yksT) d{naY 


(92) 


with the force dependent diffusion coefficient. 

One can check that the self-similar scaling solution of 
Eq. (l92|) is consistent with the average drift of the tagged 
segment in strong force regime discussed in Sec. HU As¬ 
sume that at time s, the tension gets transmitted up 
to m(s)-th segments from the pulled end. Requiring 
Eq. (1921) to be invariant under the scale transformation 
t ^ st and n —>■ n*{s)n, one obtains the dynamics of the 
tension front n*(s), which is given by Eq. (Hg. 

Note that the above stretching process can also be ana¬ 
lyzed by a different, but related nonlinear diffusion equa¬ 
tion, which describes the time evolution of the segment 
line density field [ 13 , [ 13 , . 


Appendix B 

Fluctuation-response relation in mode space 

We calculate the response function Xp(t, s) = 
S{Xp{t))/SFp{s) and the correlation function Cp{t,s) = 
{AXp{t)AXp{s)) without assuming tg —t —oo. We first 
consider the case with unchanged spring and friction con¬ 
stants, which applies to the Rouse model, and the more 
general case with SA and His in the weak force regime. 


Weak force regime 

Response function — Upon time derivative of Eq. (|44l) 
and taking ensemble average over the noise sequence 
Zpit), the response function is obtained as 

Xp(t, s) = + —5{t - s). (93) 

Ip 7p 

Correlation funetion — The time correlation of 
AXp{t) = Xp{t) — {Xp{t)) can be decomposed as 

Cp{t, s; to) = (t, s) + {t, s; to), (94) 

where the first is the stationary part invariant with re¬ 
spect to the time translation, i.e., cjf^\t,s) = Cp^\t — 
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s,0); 

s) = + —keTSit - sj95) 

7p 7p 

and the second is the excess due to the non-stationarity 
of the process; 

^ 7p 

7p 

where we add an auxiliary argument tg to indicate the ini¬ 
tial time. One can verify the FDT (HSl) . provided that the 
process is stationary, i.e., to —t —oo. Note that the excess 
part dMl) identically vanishes, when the equi-partition 
condition 

(x2(to)) = ^ (97) 

kp 

holds for each of p > 1 modes at t = to, where the aver¬ 
aging is taken over the probability distribution of Xp at 
t = to- 


For s < t < Tfo case, the response function is the same as 
that for the weak force regime (Eq. (IM1) 1. For tjq < s <t 
case, it again takes the same functional form as Eq. dMD 
with the replacement {'jp,kp) {-^p^\kp'^). Only for 
the case s < Tfo < t, the stationarity in the response 
function breaks down, and there appears an auxiliary 
argument Tfo- 

Correlation function — From the solution of Eq. (ED), 
the fluctuation in the velocity AXp{t) = Xp{t) — {Xp{t)) 
is obtained as 


Cpjit^s-,Tfo) = C^'f {t, s) + C^j\t, s;TfQ), ( 100 ) 


f fo dsxp(t, s)Z(s) - Xp(0)^e-(^»/"^»^* (r/o > t) 

AXp(t) = I fl,dsXpy^"^*\t,s;Tfo)Zy\s) 

-AXp(r/o)^e-('="V7f^>)(t-rpo) (t>rfo). 

(99) 

From this, one can calculate the correlation of the ve¬ 
locity fluctuation Cpj{t, s;Tfo) = {AXp{t)AXp{s)), and 
obtain the followings; (i) For s < t < Tfo, the weak force 
regime applies, so it is given by Eqs. (IMl) - (IMl) : (ii) For 
Tfo < s < t, it can again be decomposed as 


Strong force regime 

When the polymer with SA and/or His is stretched 
strongly, the calculation becomes slightly complicated 
due to the time dependence of spring and friction con¬ 
stants. 

Response function — From the solution of Eq. (ED: the 
response function is obtained as 


where Cy*f{t,s) and (t,s-,Tfo) take the same 

functional forms as those in the weak force regime 
(Eqs. ((551) and (1551) . respectively) with the replacement 
{Xp,kp,to) -)► {xy\ky\Tfo)-, (hi) For s < r/o < t, it 
becomes 


-i{ex) , 


Xpj{t,S-,TfQ) 

Xp{t,s) is<t<Tfo) 

.g-(fcp/7p)(i-/o-s)-(fe),^V7/^’)(‘-T/o) ("g < <t)~ 

g-(4^V7/^')(i-s) _|- '^5{t — s) {Tfo < S <t). 

7p 

(98) 


Cpj{t,s;Tfo) = -M^!^g-(fcp/7p)(r/o--)-(4^V7"b(t-^/o) 

7p7F 


■M) 


/,(/) 


7p7p^^ 


(fcsT - kp{Xp{0y)) 

( 101 ) 
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